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ABSTRACT 

We describe further observations of QSO J0906+6930, a z=5.48 blazar likely to be detected 
in 7-rays. New radio and X-ray data place significant constraints on any kpc-scale extension of 
the VLBA-dctected jet. Improved optical spectroscopy detects absorption from an intervening 
galaxy at z=1.849 and raise the possibility that this distant, bright source is lensed. We combine 
the new data into an improved SED for the blazar core and comment on the Compton keV-GeV 
\^ • flux component. 



(N 



Subject headings: galaxies: jets - quasars: general 



1. Introduction 

in 

Q0906+6930 (=GB6 J0906+6930=SRM J090630. 74+693030.8) was discovered in a survey of radio- 
\Q \ bright, flat spectrum sources chosen to be like the EGRET blazars (Sowards-Emmerd, Romani & Michelson 

2003; Sowards-Emmerd, et al. 2005). Follow-up observations confirmed a large z « 5.5 redshift and found 
evidence for a very compact pc-scale jet with the VLBA (Romani et al. 2004). The source is quite radio- loud 
with an R — /5GHz//440nm(rest) 500 — 1000, supporting the blazar interpretation. These data showed 
several peculiar properties. For example, the radio spectrum appeared to steepen above 10 GHz, while the 
jet component still appeared to be inverted. Also, the source shows a relatively large 1350A continuum 
flux and large kinematic widths for the emission lines. These suggest a large ~2x 10 9 M Q black hole mass 
[note that an error in the \F\ luminosity quoted in Romani et al. (2004) implied M > 10 10 M Q ]. While 
this source was, at best, a weak background enhancement in the EGRET survey, the prospects for detection 
with GLAST seem quite strong. This is particularly interesting as observations of a cut-off in the blazar 
spectrum above ~ 10 GeV can be used to probe absorption by light produced at the peak of star formation 
(Chen, Reyes & Ritz 2004). 

We report here on further observations of Q0906+6930, which support its identification as a high-z 
blazar, test the nature of the jet component and probe it's status as a high mass, high luminosity source. 



2. New Observations 

Several radio-loud QSO at high redshift show kpc-scale jets so our first objective was to constrain such 
arcsecond-scale emission. At the blazar redshift, 1" = 6.19kpc (for an O m = 0.24, Hq — 73km/s/Mpc flat 
cosmology, assumed throughout; Spergel et al. 2006). Perhaps the most convincing case of such jet emission 
is still that of QSO GB 1508+5714 at z=4.3, which is detected both in the X-ray (Siemiginowska et al. 2003) 
and radio (Cheung 2004) bands. Its large observed fx/ fit ~ 10 2 supports a scenario where the jet electrons 
Compton up-scatter CMB seed photons, which have a large increase in energy density, u cx (1 + z) 4 at high 
redshift (Schwartz 2002). 



-2- 



2.1. X-ray imaging 

We obtained a 30 ks Chandra observation of the blazar on 2005 July 1 positioned at the best focus of 
the ACIS-S3 backside-illuminated chip. The observations were standard full-frame 3.2s TE exposures in 
very faint (VF) mode. No strong background flaring was seen during this observation so that all of the data 
could be included in the analysis. The data were processed using standard algorithms in CIAO version 3.2 
and CALDB version 3.0.0. The blazar was clearly detected, providing ~ 435 counts, at a position coincident 
with the radio core. With this modest flux, pile-up is under 1%, and so we can ignore it in the subsequent 
analysis. The X-ray image is consistent with a PSF computed for the blazar spectrum; the largest excess 
(at ~ 1.5", position angle ~ 115°) contains 4 photons and is not significant. Thus, on ~ 1 — 10 arcsecond 
scales, any extended (jet) emission must have < 1% of the core flux. A number of coronal emission field 
stars and some unidentified objects (presumably background AGN) were detected in this integration, but 
none showed notable peculiarities, so we do not discuss these further. 
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Fig. 1. — CXO ACIS-S3 spectrum of Q0906+6930, with the best-fit absorbed power-law model. The second panel 
shows the 0.5-8 keV CXO image. 

The X-ray spectrum (figure 1) is adequately described by a simple absorbed power-law, with a low 
absorption N H = 5.4± 2 x 10 20 cm -2 , a photon index T = 1 .6±0.1 and an unabsorbed flux f x (0.2 — 10keV)) = 
1.47 ± 0.25 x 10~ 13 crgcm~ 2 s _1 . Although the data show a small excess to the best fit power-law model 
around 1 keV, an attempt to fit a red-shifted Fc line produced no significant detection. The upper limit to 
the line flux is ~ 1.5 x 10~ 15 ergcm~ 2 s _1 , or < 1% of the continuum flux. These spectral parameters are 
quite typical of lower redshift blazars. The measured absorption is consistent with the Galactic column in 
this direction, ~ 4.4 x 10 20 cm~ 2 from the COLDEN tool. We can also compare the X-ray properties with 
those of a set of radio- loud high- redshift QSO recently reported by Lopez et al. (2006). Our observed T 
and inferred aox [250 nm/2keV] (rest) w — f.3, arc in excellent agreement with the objects measured in this 
paper. We conclude that this is a fairly typical radio-bright quasar at high redshift. 
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2.2. Radio imaging 

To search for a larger scale counterpart of the VLBI jet reported in Romani et al. (2004), we also 
obtained VLA images with the A array at 1.4- , 4.8- and 8.4 GHz on 2004, October 7. The on-source times 
were 10m, 30min and 50min respectively. These data were calibrated and mapped using standard AIPS 
routines and were found to be adequately described by a simple point source model. The final map noise 
was within a factor of <~ 2 of the thermal noise in each case. Measuring the rms background in an annulus 
between the beam size and ~ 5", we obtain limits on any jet component (Table 1). The largest map fluxes 
in the searched regions are consistent with statistical fluctuations for these rms vales, so we can take upper 
limits on the jet flux along the position angle of the mas jet to be ~ 2x the local rms, while the limits on any 
extended component are about two times larger. For the X-band image this search extends within ~ 0.4" of 
the blazar core, while for Q-band in the B-array, we searched down to 0.1". 

Since no extended emission is seen in the A-array data, we used follow-on B-array observations on 2005 
March 04 to obtain a six frequency, single epoch radio spectrum. The on-source total scan times are listed 
in Table 1. The data were calibrated by observations of 0542+498 (=3C 147) at each frequency. For all 
measurements, the statistical error in the flux density estimate was small. We add this in quadrature to 
the flux scale uncertainties described in the VLA calibrator manual (2% at L, C and X and 5% at U, K 
and Q) to obtain the flux density errors in table 1. For K and Q residual pointing errors and uncorrected 
atmospheric absorption may contribute additional error. In all cases a single unresolved component fit 
provided a satisfactory description of the Q0906+6930 data. In Table 1, we also give RMS fluctuation levels 
in an annulus about the point source. The local rms limits to extended flux for the three lower frequency 
bands are from the deeper, higher-resolution A-array images. 

Table 1: Radio Core Fluxes and Jet Bounds 
Band Freq(GHz) On-Source time (min) Core Flux (mjy) Error (mjy) Jet Annulus RMS (mjy) 



L 1.44 8 92 2.0 O.IO 1 

C 4.88 8 114 2.4 0.05 f 

X 8.44 12 136 3.1 0.13^ 

U 14.96 24 129 6.6 0.23 

K 22.49 15 83 4.2 0.22 

Q 43.33 30 43 2.2 0.26 
RMS from 2004/10/04 A-array maps 



2.3. Optical Spectroscopy and Imaging 

To further investigate the absorbed optical spectrum reported in Romani et al. (2004), we obtained 
higher resolution spectra using the Hobby»Eberly telescope (HET; Ramsey et al. 1998) Marcario Low Reso- 
lution Spectrograph (LRS; Hill et al. 1998) using the G3 VPH grating, which provides 1.88A/pixel dispersion 
over the range AA6250A - 9100A. The integrations used a 1.5" slit, giving an effective resolution of 5.5A. 
Standard queue observations were made on 2005 February 9 and 14 and 2005 May 9-10; conditions were 
generally sub-par with poor transparency or variable seeing. After standard IRAF processing and combi- 
nation using the relative S/N, seven exposures totaling 8600s of integration were assembled to produce the 
spectrum in figure 2. Here the new G3 spectrum is plotted with a heavy line over the lower resolution Gl 
spectrum of Romani et al. (2004). The excellent correspondence between the fine structure in these spectra, 
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especially in the Lya forest region, attests to the stable calibration. The combined spectrum has a S/N per 
resolution element of ~ 16 in the unabsorbed continuum. 

These higher resolution data allow us to check the reality of absorption-line systems suggested in the Gl 
spectrum. Two metal line systems can be seen in the inset of figure 2, which shows the spectrum redward of 
Lya. In particular, we confirm a likely C IV doublet at z = 4.698, with a rest equivalent width (REW) for 
each component of Wq = 0.6 ± 0.13A. This system has an associated, apparently damped, Lya absorption 
and a minimum at Ly/3. More interesting is an apparent strong Mg II absorption at z = 1.849 which was 
confused by possible N V emission in the lower resolution data. The REW of the 2795A line is 1.9 ± 0.07A 
and the line spacing and doublet ratio R = 1.2 support this identification. In our restricted spectral range 
no strong confirming line is expected. 

The presence of intervening absorption line systems raises the possibility that the blazar is magnified 
by gravitational lensing, as discussed below. To constrain this hypothesis, we examined the limited imaging 
data. During HET queue spectroscopy, 10-30s ~ white light set-up images are generally recorded through 
a 515nm long-pass filter. Six of these images had usable (< 1.7") image quality; these were S/N weighted 
and combined to produce a final image with 1.6" FWHM. This was used to get an approximate limit on any 
extended source within 5" of Q0907+6930 of R ~ 23.1. Imaging reconnaissance was also obtained using the 
MDM 2.4m RETROCAM on 2006 February 4, courtesy of J. Marshall and D. DePoy. Here 6 x 300s r images 
(~ 1.3" FHWM) were combined. Although Lya forest absorption should suppress the quasar relative to 
the foreground galaxy in the r band, in the combined image the blazar was again stellar, with an apparent 
magnitude of r — 21.9. We derive an upper limit of r > 23.3 on any nearby galactic companion. 

3. Spectral Energy Distribution 

We have assembled these new flux measurements into a (non-instantaneous) spectral energy distribution 
(SED), figure 3. The plot shows the 04/03/05 simultaneous radio fluxes. We do not see any strong evidence 
of radio variability, but it should be noted that the X-ray flux (shown with a dashed lines for the continuum 
uncertainty range) was observed ~ lOOd later. The optical spectroscopy measurements span several years, 
although again we do not see variability. With a rather small range of unabsorbed continuum available, the 
optical spectral index is quite uncertain (dashed lines, extended a a decade for visibility). This situation 
would be greatly helped by near-IR spectroscopy. Finally, the EGRET upper limit is, of course, from 10-15 
years prior to the current observations, and blazars are known to be highly variable in the 7-ray band. The 
low significance detection during one viewing period (Romani et al. 2004) should be seen as consistent with 
the mission- averaged upper limit. 

As in Romani et al. (2004) we compare the new SED with that of the well-studied EGRET blazar 
3C 279 [inset: SED from period P5b (Hartman et al. 2001), red-shifted to z=5.48], and with synchro- 
Compton models. The Q0906+6930 spectrum is sparse, but notably shows a significantly (> 3x) more 
luminous optical spectrum. In contrast to 3C 279 the optical emission shows high equivalent-width lines 
from the broad-line region (BLR), and is thus not dominated by synchrotron continuum. For the Compton 
components of the SED, our new X-ray detection confirms a rising vF v spectrum, likely up-scattered by the 
synchrotron-emitting electrons. This is somewhat (~ 50%) fainter than the 3C 279 X-ray luminosity, which 
is reasonable in light of the lower radio synchrotron luminosity. However, the possible GeV detection would 
require an external Compton component ~ 30 x brighter than that of 3C 279. 

By computing representative synchro-Compton models, we see that the brighter disk emission and a 
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luminous 7-ray Compton component may be related. In the inset we show a model for 3C 279 computed 
from a code kindly supplied by M. Bottcher. This code computes synchrotron and Compton emission in an 
evolving cylindrical jet outflow. It also models a host accretion disk and computes the Comptonized flux 
from these soft disk photons, both directly from the disk and disk photons scattered from the broad line 
region. The code is described in Bottcher, et al. (1997) and Bottcher & Bloom (2000). To generate the 3C 
279 model, we use explicitly the parameters listed in Hartman et al. (2001) for this epoch. The 'jet' is a 
6 x 10 16 cm radius blob injected 0.025 pc from the 1.5 x 10 8 M Q hole, moving cylindrically at 2° to the line of 
sight with bulk T = 13. The jet contains 35cm -3 each of e ± , with a power-law energy distribution of index 
p = 3.0, extending from 71 = 6 x 10 2 to 72 = 1 x 10 5 and 1.5 G of magnetic field. The hole is surrounded by 
an accretion disk of luminosity 10 erg/s and a broad-line region extending from 0.1-0.4pc. To reproduce 
the model shown in (Hartman et al. 2001), it was necessary to increase the BLR Thomson depth to t = 0.06. 
These parameters gave the total SED curve shown in the inset. 

We know that Q0906+6930 has a higher mass hole (~ 2 x 10 9 M Q ) and a brighter accretion disk 
(Ld > 3 x 10 46 crg/s). We made these two changes to the parameters and re-computed the SED, as shown 
in the main part of Figure 3. Under the full SED we show the individual components: (from low to 
high frequency) the jet synchrotron emission, the disk thermal emission, the jet self-Compton emission, the 
external Compton emission from the disk and the Comptonized emission from disk photons scattered off the 
broad-line region. The spectra were computed at z = 0.538 and then the 3C 279 data and the two models 
were red-shifted to z=5.48 using concordance flat cosmology parameters. These models illustrate that with 
a higher mass hole and a brighter, cooler disk, brighter external Compton (7-ray) emission can be expected. 
Of course, such models have many parameters and fits to more densely sampled SEDs in a variety of flare 
states are needed for any predictive power. 

A standard definition of radio-loudness is the ratio R = /sGHz//440nm(rest) (Kellerman ct al. 1989). 
Although we do not have direct observations at 770 MHz or 2.9/x, we can extrapolate the observed optical 
and radio spectra to estimate R « 540 — 970 where the range is due to the rather uncertain optical continuum 
spectrum. 

We should note that both the 3C 279 and the Q0906+6930 models safely under-predict the radio 
spectrum. This is because the model simulates only the base of the jet emission, omitting larger scales where 
cooling, expansion and injection of additional populations produce the cm- wavelength radio components. 
In particular, the local maximum in the Q0906+6930 radio spectrum at ~ 15GHz (~ 100GHz rest-frame) 
suggests a distinct radio jet component. It should be noted that the 43 GHz flux measured here (on kpc 
scales) is in good agreement with the VLBA core measurement. Thus at most 10% percent of the radio flux 
lies in the ~ 5 — 10 pc resolved VLBA jet and the cm- wavelength radio component is produced at sub-pc 
scales in the unresolved VLBA core. It would be useful to check if the SED rises toward a synchrotron peak 
in the far-IR/sub-mm region as expected. IR photometry of the blazar will also be helpful in constraining 
the level of the synchrotron SED and confirming the dominance of disk emission in the optical. 

4. Conclusions 

The new data generally confirm the SED picture described in Romani et al. (2004) and show that at least 
the X-ray Compton component is present in this blazar. The lack of any resolved kpc-scale jet component 
suggests that all emission observed is from the core (or a compact pc-scale jet). Also, this does not support 
the Schwartz (2002) picture where large-scale X-ray jets are produced by up-scatter of the CMB. In this 
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model, the oc (1 + z) A increase of the CMB energy density counters the decreasing jet surface brightness; at 
z > 4 jet emission is thus expected to exceed core emission. With < 1% of the core X-ray flux in a resolvable 
jets, our results extend those of Lopez et al. (2006), suggesting core photons, rather than CMB photons 
dominate the flux up-scattered as Compton X-ray emission. 

One peculiarity noted in Romani et al. (2004) was the relatively large optical and radio flux; we noted 
the possibility of lensing of the QSO core. This should be re-visited, since the strong intervening metal line 
absorbers suggest galactic-scale masses along the line-of-sight. The most important is the Mgll system at 
z=1.849, which at Wq = 1.9 A., likely represents a classical ~ spiral galaxy (Churchill et al. 2005) which 
should be a Lyman Limit System (LSS) absorber, within an impact parameter 

R = 38h- 1 (L K /L* K ) 015 kpc, 

for an L» galaxy at the absorber redshift, i.e. within 6.2" for our assumed cosmology. Our limited direct 
imaging does not detect such a galaxy, but the results are not very constraining. At this z, the cosmic age 
is ~ 3.6 Gyr and an L* Sa formed at z=5 in our cosmology would have r ~ 23.9 for a Bruzual & Chariot 
(1993) model, i.e. about 2x fainter than our imaging r ~ 23.3 limit. Sbc types would be even fainter. 

Can such a galaxy lens our blazar? The expected image angular scale for a lens mass M is 

e = [AD LS GM/(D s D L c 2 )] 1/2 

where one uses angular diameter distances and D LS = D s — (1 + z L )D L /(l + z s ). The geometry is favorable 
for lensing and the lens scale is then 1.35 x 10 _6 M 1 / 2 arcsec with M in M Q . We see no evidence for a double 
image of the blazar on arcsecond scales, either in the optical or the radio, so this precludes classical macro 
lensing by an L* galaxy. Indeed, the core is unresolved in our VLBI maps, where the half-power beam width 
at 43GHz was 0.55x0. 30 mas, so we can place an upper limit on an effective lens mass of M < 1.4 x 10 5 M Q . 
The lack of obvious variability on year time-scales suggests a not very constraining lower mass limit of 
M > lO _3 w| oo M0 for a lens galaxy velocity dispersion of ~ 300f3ookm/s. We may infer a slightly tighter 
limit if we note that the large BLR equivalent widths imply that the broad line region must be lcnscd along 
with the core. Naively applying the continuum BLR radius correlation of Kaspi et al. (2005) to our optical 
continuum luminosity ^aI 5100 ^ ~ 3 x 10 46 erg/s we conclude 

/ \0.69 

Rblr « 5.8 x 10 16 (AF A | 5ioo ^/10 44 erg/sJ » 2.9 x 10 18 cm. 

This in turn implies a minimum lens mass M > 10 _2 Mq. We infer that macro-lensing does not amplify 
Q0906+6930, but that micro- or milli-lensing are possible. On the other hand, it should be noted that the 
probability of an unrelated intervening galaxy is not small. Since dN/dz(W > 1.9 A) w 0.13 over the redshift 
range covered by our G3 spectrum (Nestor et al. 2005), we expect 0.06 such absorbers. Detection of one 
Mgll systems is thus not very improbable and so we cannot infer from this that lensing amplification has 
brightened Q0906+6930. 

We conclude by re-iterating that our new SED measurements support the picture of Q0906+6930 as a 
blazar with a bright 7-ray Compton component emitting at z=5.48. This Compton radiation must traverse 
the peak of star formation at z ~ 2 — 3, where optical/UV emission can attenuate the >10GeV photons, 
just as the IR background can attenuate the TeV emission of blazars at lower z, e.g. Dwek, Krennrich 
& Arendt (2005) and references therein. Of course, there is also attenuation from the host photon field 
and, in principle, from the intervening z=1.89 galaxy. These contributions will make it difficult to extract 
constraints on the extragalactic background photon field from this (or any one) object. However, statistical 
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studies of high red-shift blazars with GLAST (Chen, Reyes & Ritz 2004), should still be able to extract 
global constraints on the extragalactic background light and its evolution. 
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Fig. 3. — Multiwavelength spectral energy distribution for Q0906+6930, from non-simultaneous data (open circles). 
The optical continuum flux is plotted with a line (short dash) extended to show the slope. The X-ray spectrum 
(solid line) shows our CXO continuum fit. For these two components dashed lines show the range of spectral index 
uncertainty. For the 7-ray band we show the brightest epoch EGRET flux, the EGRET mission-averaged upper 
limit and the GLAST sensitivity of the first year sky survey. A synchro-Compton model is plotted (see text). The 
individual component curves (synchrotron emission, synchrotron self-Compton emission, disk emission [dashed line] 
and Comptonized emission directly from the disk and after scattering in the broad line region) are shown. For 
comparison, the inset shows the SED of 3C 279 during period P5b Hartman et al. (2001), shifted to z=5.48, along 
with its synchro-Compton model (see text). 



